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ABSTRACT 

In a practical design of an electric -propulsion 
space vehicle, it is likely that some components, 
such as a solar -cell array, will be exposed in some 
degree to neutral atoms and charge -exchange ions 
leaving the thrustors . The purpose of this paper 
is to describe the energy and angular distributions 
of such particles emerging from a 15 cm diameter 
mercury electron-bombardment thrustor. An analyti- 
cal approach was used to determine the particle ef- 
fluxes. The results indicate that the erosion and/ 
or coatings due to charge -exchange ions or neutral 
particle effluxes would be tolerable for several 
years if components were located outside the pri- 
mary ion beam and at distances of 75 cm or greater 
from the thrustor. The value of 75 cm was chosen 
to simplify the analytical solution. At distances 
much closer than 75 cm a more detailed study would 
be required because the thrustor can no longer be 
accurately represented as a point source. Approxi- 
mate scaling relations are given that should allow 
rough determination of particle effluxes for vari- 
ous combinations of nominal operating parameters 
other than those used in this study. 


INTRODUCTION 

Electric propulsion is becoming more attrac- 
tive for a variety of spaceflight missions as a re- 
sult of improvements in thrustor performance and 
component technology. For example, an interesting 
design study has been made recently to investigate 
the feasibility of electric propulsion systems for 
unmanned probes (ref. l) . One of the possible 
problems found in the course of such design studies 
concerns the relative placement of thrustors with 
respect to the solar -cell arrays or other space- 
craft components (ref. 2). If the spacecraft is 
designed for maximum compactness and minimum weight, 
it is likely that the solar-cell array will be ex- 
posed to some amount of neutral atoms and/or ions 
leaving the thrustor. The major directed part of 
the ion beam is easily avoided. However, the "edge" 
of the beam is not well defined due to the presence 
of charge -exchange ions in the thrustor exhaust. 

In recent long duration tests, neutralizers were 
mounted in regions outside the major directed ion 
beam and they experienced appreciable ion erosion 
(ref. 3). The extent of the damage was found to 
decrease as the radial distance from the beam axis 
increased. Thus, the problem arises of defining 
regions in which possible neutral particle coatings 
and/or ion erosion will be at an acceptable level. 
Of course, acceptable levels must be determined 
from considering the mission duration, the spectrum 
of particles leaving the thrustors, and the proper- 
ties of the spacecraft structure hit by these par- 
ticles. 

The purpose of this paper is to describe the 
theoretical energy and angular distribution of neu- 
tral atoms and charge -exchange ions leaving a 
15 centimeter diameter, mercury electron bombard- 




ment thrustor. Because the magnitudes of these 
quantities are small, accurate experimental measure- 
ments are very difficult. However, the processes 
giving rise to such effluxes are fairly well under- 
stood so that an analytical approach may be used. 
Although the results obtained are not precise, they 
are adequate for initial spacecraft design calcula- 
tions. In addition, they serve to indicate which 
aspects of the problem may deserve more detailed 
study. 


Model and General Approach 

The thrustor used as the basis of this study is 
described in reference 3 and shown in figure 1. At 
normal operating conditions the beam current from 
this thrustor is 250 mA, extracted with a screen 
voltage of 3 kV and an accelerator voltage of -2 kV. 
With 511 holes in the electrodes, the average cur- 
rent per hole is about 0.49 mA. The extracted cur- 
rent distribution across the thrustor is such that 
at the center the maximum current per hole is about 
1.0 mA, or about twice the average value. The cur- 
rent extracted from the outer perimeter holes is 
about 0.25 mA per hole. The thrustor is normally 
operated at a propellant utilization efficiency of 
80 percent. Anode wall temperature is about 500° K, 
resulting in a neutral atom temperature in the dis- 
charge chamber of about the same value. 

For the analysis, the possible charge -excharge 
ions were grouped as follows: 

(1) those leaving the thrustor with dispersion 
angles between 0° and 20° (considered to be 
within the primary ion beam) 

(2) those leaving the thrustor with dispersion 
angles between 20° and 90° 

(3) those returning to the thrustor and imping- 
ing on the accelerator electrode 

(4) those produced downstream of the neutrali- 
zation plane and not returning to the thrustor 

The ions in groups ( 2 ) and (4) were of prime 
interest for this study, since these ions may im- 
pinge on surfaces downstream of the thrustor. Fig- 
ure 2 shows typical trajectories of the four groups 
of charge -exchange ions and the approximate region 
of origin. The zero volt equipotential line, which 
falls between the screen and accelerator, was re- 
garded as the downstream boundary of the formation 
region for the group (2) ions. This was because any 
ion formed at a negative potential would return to 
the thrustor. 


The electrode configuration shown in figure 2 
represents one hole of the thrustor extraction sys- 
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that point. The angle between the centerline and a 
tangent to the trajectory at the point where the ion 
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enters the neutralized exhaust plasma was used to 
determine the angular distribution of the charge - 
exchange ions. The formation of the group 4 
charge -exchange ions in the region downstream of 
the neutralization plane was treated as a separate 
problem and will be discussed later. 

The charge -exchange analysis was made using an 
IBM 7094 digital computer program (ref. 4) and an 
electrolytic tank analog (ref. 5) in a combined 
manner described in detail in reference 6. Briefly; 
the technique involves determining the Poisson po- 
tential distribution of the desired model using the 
digital computer program, then establishing this 
distribution in the electrolytic tank. (With this 
method the potential distribution of an axially 
symmetric model is readily established in a simple 
two-dimensional (flat bottom) tray of the electro- 
lytic tank analog.) Once the proper distribution 
has been established, the exhaust beam is arbitrar- 
ily divided into a number of mesh squares- Charge 
exchange ion trajectories are then started with 
zero initial velocity from the centers of these 
mesh squares and traced on an X-Y plotter. The lo- 
cal charge -exchange ion production rate is then 
calculated as follows. The charge -exchange ion 
production rate in an element of volume is 

N » Qnn o AL a (l) 

where Q is the cross section for charge -exchange 
(the value 6X10“-^-^ cm^ was used for mercury) . (All 
symbols are defined in appendix A.) The elements 
of volume in this case are toroidal with square 
cross section (a * 2rtr AL, where r is the average 
radius of the toroid) . 

The primary ion arrival rate n, was calcu- 
lated by a technique similar to the one used in de- 
termining charge contribution of flowing ions to 
elements of volume in reference 5. In brief, the 
primary ion flow is divided into stream tubes (the 
ion trajectories the stream tube boundaries) . 

The current flowing in each stream tube is known 
from the digital computer solution. The ion ar- 
rival rate at a particular volume element is thus 
determined by summing the contributions from all 
the stream tubes intercepted by the volume element. 

The neutral atom density, riQ, was assumed to 
be uniform in the interelectrode region and extend- 
ing downstream as far as the neutralization plane. 

A value was obtained from the relation 

_ H 0 ta H 0 ta _ Ho 8 - , ^ 

^0 “ V = A AL “ Av 

where is the rate at which neutrals leave the 

screen aperture, A is the area of a screen hole, 
and v is the r.m.s. velocity of a neutral atom 
at 500° K. The assumption of a uniform neutral 
density for the entire extraction system implies 
that the propellant utilization efficiency per hole 
varies since the current density per hole is not 
uniform. The total charge -exchange ion production 
rate in the interelectrode region was determined by 
summing all the local contributions calculated by 
equation (l) . Results were obtained for a few typ- 
ical cases and are presented in the RESULTS AND 
DISCUSSION section. In order to generalize the re- 
sults to cover other cases, the effect of changing 
various thrustor parameters on the charge -exchange 
ion production was examined. The resulting scaling 
relations are given in appendix B. 


RESULTS AND DISCUSSION 
Charge -Exchange Ions 

The charge -exchange ion analysis was performed 
for three different values of current extracted 
from a single hole in the electrode. The current 
values were approximately chosen to represent ex- 
traction from a typical outer hole, midsection hole, 
and central hole of a thrustor operating at a beam 
current of 250 mA. The specific values used for 
individual hole currents were 0.230 mA, 0.435 mA, 
and 1.034 mA, respectively, with corresponding pro- 
pellant utilization efficiencies of approximately 
65, 80, and 90 percent. Figure 3 shows the results 
obtained for group 2 charge exchange ions for the 
high current, central hole case. Part (a) of the 
figure shows some typical charge -exchange ion tra- 
jectories. In part (b) the charge -exchange ion 
formation region is shown divided into subregions 
which produce ions in an indicated range of disper- 
sion angles. No subregion is shown producing ions 
with dispersion angles above 44°. This lack of de- 
tail is a result of the mesh size used. Actually, 
there will be a distribution of ions between 45° 
and 90° dispersion angles that originate in the 
downstream half of the two subregions near the axis 
on the right. These ions, however, would have very 
low energies since the zero volt equipotential is 
very close to the region at that point. Also not 
indicated, due to the coarseness of the mesh used, 
is a narrow region near the centerline of the hole 
where the ions formed will have dispersion angles 
less than 20° (group (l) type). In part (c) of fig- 
ure 3 the charge -exchange ion formation region is 
divided into subregions that produce ions in a par- 
ticular energy range. 

Results for the midsection and outer hole were 
somewhat similar in appearance to figure 3. It was 
interesting to learn, however, that the maximum 
dispersion angle of the primary ion beam was largest 
for the outer hole. This is a result of "crossover" 
trajectories (ref. 5) rather than simple beam 
spreading. It was previously thought that the pri- 
mary ion beam boundary existed at approximately a 
15° dispersion angle; however, this study and the 
long duration tests on thrustor s indicate that a 
low primary ion flux continues out to about a 20 u 
dispersion angle. 

The results of the analysis are presented 
graphically in figures 4 and 5. Figure 4 shows the 
average number of group (2) ions within a particular 
range of dispersion angles. In all three cases most 
of the ions leave the thrustor at relatively low 
dispersion angles. Figure 5 shows group (2) ion 
energies versus dispersion angle. The solid curves 
represent the maximum energies determined and the 
dashed curves are an extrapolation of the solid 
curves to the 90° dispersion angle. The shaded re- 
gions cover the areas in which all the values of 
this analysis occurred. Charge -exchange ions will 
be produced with energies filling the remainder of 
the area under the curves, but they will be rela- 
tively small in number. 

The groups (l) and (3) charge -exchange ions 
were determined in number only (i.e., specific dis- 
persion angles, or, impact points on the accelerator 
for the group (3) ions were not recorded) . The to- 
tal number of these charge -exchange ions was calcu- 
lated. using equation (l) . The results obtained for 
the first three groups of charge -exchange ions are 
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summarized in table I. From the table it can be 
seen that the total charge -exchange ion current 
from groups (l), (2), and (3) ions was about 0.4 
percent of the total beam current. Of this total 
charge -exchange ion current, 84 percent impinged 
on the accelerator (group (3)), 10 percent escaped 
with the primary beam (group (l)), and 6 percent 
escaped with a dispersion angle between 20° and 90° 
(group (2)). Figure 6 shows the total charge- 
exchange ion flux for group (2) ions as a function 
of dispersion angle. In preparing this graph it 
was assumed that the thrustor could be regarded as 
a point source. A distance equal to about five 
thrustor diameters (75 cm) seems sufficient for 
this purpose. 

To calculate the group (4) charge -exchange 
ions (those produced in the region downstream of 
the neutralization plane) , the following assump- 
tions were made: (l) neutrals are emitted uni- 

formly over the thrustor diameter and follow a co- 
sine law of divergence, (2) the ion flux leaving 
the accelerator grid is paraxial, remains paraxial, 
and is of uniform density. The simplifying assump- 
tion of no beam spreading results in negligible 
error of total number of char ge -exchange ions pro- 
duced because: (a) neutral density does not change 

appreciably in the region between 0° and 20° beam 
dispersion, and (b) the same number of primary ions 
are considered in either case. 

A radial neutral density distribution was cal- 
culated for several distances downstream of the 
thrustor using equation (6) of reference 8. The 
result of these calculations are shown in figure 7. 
The densities have been normalized to the value of 
neutral density at the neutralization plane. The 
neutral flux gradients that would occur near the 
thrustor because of the grid structure were assumed 
to blend into a smooth profile approximately 3 cen- 
timeters downstream. 

Once the neutral density distribution is de- 
termined, the calculation of the charge -exchange 
ion formation rate at various locations can be 
made using equation (l) . The calculation was made 
using the thrustor operating conditions of 80 per- 
cent propellant utilization efficiency, 0.250 am- 
pere beam current, and 3000 volts net accelerating 
voltage. A value of 6X10~15 cm 2 was used for Q. 
The result is plotted in figure 8 as the number of 
ions per cm 2 per second passing through the wall of 
a cylinder equivalent in shape to the assumed pri- 
mary ion beam (see sketch of cylinder on fig. 8). 
Each charge exchange ion was assumed to travel 
radially and normal to the axis at the beam. This 
assumption seems valid for the reasons given below. 
The charge -exchange ion has only random thermal 
(~500° K) energy when it is formed. Any electric 
fields existing in the ion beam will dictate the 
direction of travel of the charge -exchange ions. 
Reference 9 indicates that the predominate (al- 
though small) electric field in the ion beam will 
be radial. The energy gained by a charge -exchange 
ion leaving a well neutralized beam is probably 
less than about 50 volts. Figure 8 indicates that 
the maximum number of group 4 ions will be found 
near the thrustor. The decreasing number found 
further away from the thrustor is, of course, due 
to the diminishing neutral efflux. 


Effects on Spacecraft Components 

In order to give an idea of the possible af- 
fects of the neutral atom effluxes and group (2) 
and (4) charge -exchange ions on spacecraft compo- 
nents two different shaped surfaces located down- 
stream of the thrustor were considered as intercept- 
ing these effluxes. The first surface was a cylin- 
der, concentric with the primary beam, and with a 
diameter of 150 centimeters or 10 times the thrus- 
tor diameter (see fig. 9(a)). The primary ion beam 
was assumed to spread with a dispersion angle of 
20°. The length of the cylinder was chosen so as 
to just intercept the primary ion beam. The second 
surface was a plane, normal to the thrustor axis, 

75 cm downstream and with a hole to allow the pri- 
mary ion beam to pass through (see fig. 9(b)). 

A distance of 75 centimeters seemed adequate 
for this purpose so that the thrustor could be re- 
garded as a point source. Figure 10 shows the dis- 
tribution of the total ion flux of group (2) and 
group (4) charge -ex change ions arriving at the cy- 
lindrical surface as a function of dispersion angle. 
The highest flux occurs at low dispersion angles 
going down to a minimum flux between 45° and 50° 
dispersion, and then rising again at high disper- 
sion angles due to the increasing contribution of 
group (4) ion fluxes. 

The flux arriving at the planar surface normal 
to the thrustor axis is entirely due to the group 
(2) charge -exchange ions. Figure 10 also shows the 
distribution of this flux as a function of disper- 
sion angle. The highest flux again occurs at low 
dispersion angles, dropping to zero at the higher 
dispersion angles. 


The charge -exchange ion flux arriving at the 
surfaces used in the previous calculation can vary 
from zero to about lO^- 1 particles per cm 2 per sec- 
ond. Their energy range is from 0 to about 2.5 kV. 
In comparison, cosmic radiation flux in deep space 
is about 2.0 particles per cm 2 pex- second (refs. 10, 
11, and 12). This flux consists mainly of protons 
with energies from 10 to 10 12 MeV. On the other 
hand the proton flux in the Van Allen belts ranges 
up to 10 9 protons per cm 2 per second with an energy 
range from 0.1 keV to 700 MeV (refs. 10, 11, and 
12) . Damage to solar cells by the natural radiation 
flux is caused primarily by the penetration of the 
high energy particles into the cell itself which 
can cause a degradation in performance. The flux 
from the thrustor does not possess sufficient energy 
to penetrate a shield such as a 0.1 millimeter 
thickness of quartz commonly used to protect the 
cell from natural radiation. Rather it seems likely 
that any damage done by the thrustor efflux would 
be a result of sputtering erosion. If a sputtering 
yield of one particle per incident ion for mercury 
on quartz (refs. 13 and 14) and a flux of lCp-^ ions 
per cm 2 per second is assumed, the resulting wear 
rate is 3.8X10"-^ centimeters per second. The wear 
rate of silver, assuming a sputtering yield of 
seven atoms per incident ion (refs. 13 and 14) and 
a flux of id 1 i ons per C m 2 p er secon d is 
5.87X10“13 cm/sec. Most materials have wear rates 
between these two values. At these wear rates it 
would take 1.8X10 6 hours to wear away one-quarter 
of a 0.1 millimeter thick quartz shield. Of course, 
this wear rate would increase for components lo- 


The neutral flux was calculated by assuming a 
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uniform flux across the surface of the accelerator 
grid, a cosine distribution for the dispersion, and 
a total flow of 3.9X10 17 neutrals per second which 
corresponds to the operating conditions used in 
this study. Figure 11 shows the distribution of 
neutral particles intercepted by the cylindrical 
and planar surfaces shown in figure 9. The energy 
of the neutrals is below the sputtering threshold 
energy, therefore the most probable effect of these 
neutrals would be film deposition. Thus, it be- 
comes important to compare the arrival rate of neu- ' 
trals with their evaporation rate at the assumed 
conditions. The evaporation rate relation is 

° (2«mkT) 1//2 ^ ^ 

Assuming a temperature of 330° K for the com- 
ponent (ref. 15), this corresponds to an equilib- 
rium vapor pressure for mercury of 2X10" 2 torr. At 
these conditions, the evaporation rate is 2.73X10^ 
particles /cm 2 per second. Comparing this value 
with the values on figure 11, it can be seen that 
the possibility of forming a coating on these sur- 
faces, other than a chemisorbed monolayer, is in- 
deed remote. For cesium propellant the evaporation 
rate would be approximately three orders of magni- 
tude less. For the arrival fluxes calculated, this 
lower evaporation rate is still enough to preclude 
coating. 


CONCLUSIONS 

The results of this study indicate that the 
effects of charge -exchange ions and neutral parti- 
cle effluxes from a 15 centimeter diameter mercury, 
electron-bombardment ion thrustor on spacecraft 
components will be minimum outside the primary ion 
beam and at distances of at least 75 centimeters 
from the thrustor. Under the assumed conditions, 
silver, for example, would wear at a rate of 
3 . 8 X 10 “^ centimeters per second and quartz would 
wear at a rate of 5.d7xI0”-*-~ centime i«xt> per sec- 
ond. For most components this wear rate should be 
tolerable for missions of less than 100,000 hours. 
At distances much closer than 75 centimeters a 
more detailed study would be required because at 
these closer distances the thrustor cannot be ac- 
curately represented as a point source. 

Coating of components with neutral mercury 
atoms should not be a problem at distances of at 
least 75 centimeters since the particle effluxes 
are about seven orders of magnitude lower than the 
evaporation rate at an assumed component tempera- 
ture of 350° K. 

The approximate scaling relations given in ap 
appendix B should allow rough determination of par- 
ticle effluxes for various combinations of operat- 
ing parameters other than those used in this study. 


appendix: a 

SYMBOLS 

A area of screen hole 

a cross sectional area of volume element 

1- T> _-| J ■ . 

xv -uw. x. O/iiiiai nt v^uuo o 


AL incremental length 

m mass 

N charge -exchange ion formation rate 
n ion arrival rate 

Hq neutral atom density 

n o neutral atom evaporation rate 

P pressure 

Q charge -exchange cross-section 
r average radius of torroid 
T temperature 

t time 

V volume 

v velocity 

propellant utilization efficiency 
fig arrival rate of neutrals 

<p electric potential 


APPENDIX B 
SCALING RELATIONS 

Scaling relations were developed by consider- 
ing "equivalent cases." Two cases were considered 
equivalent if their potential distributions as de- 
termined from digital computer analysis were simi- 
lar in shape and if the location of the zero volt 
equipotential line was about the same in each case. 
The concept of equivalent cases is useful iu that 
all trajectories, both those of primary ions and 
charge -exchange ions, follow identical paths. To 
begin, several cases of interest are first run on 
the digital computer to obtain potential distribu- 
tions. Charge -exchange currents are then deter- 
mined. Now, suppose a solution is desired for a 
case that has a different combination of operating 
parameters than any of the ones previously run. 

This new case is first run on the digital to obtain 
a potential distribution. The previous cases are 
examined for one which is "equivalent" to the new 
case. Then by using the scaling relations given 
below it is possible to determine approximate values 
of the charge -exchange currents. The scaling pa- 
rameters developed in general apply only to elec- 
tron bombardment thrustors. 


Propellant Utilization Efficiency 

An expression relating the primary ion flow 
rate n, the neutral atom arrival rate Pq, and the 
propellant utilization efficiency r\ is as fol- 
lows : 


U 


0 = 


(1 - T, u ) 


n 


Now, from equations (l), (2), and (Bl) 


(Bl) 
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TABLE I. - GROUP (l), (2), AND (3) CHARGE-EXCHANGE ION PRODUCTION. 


Exhaust 

beam 

region 

Current 

per 

hole, 

mA/hole 

Ifumber 
of holes 
allocated 

Primary 

beam 

current, 

mA 

Propellant 

utilization 

efficiency 

Charge -exchange ion current, 

mA 

Group (l) 

Group (2) 

Group (3) 

Total 

Central 

1.034 

85 

87.9 

Q.90 

0.016 

0.018 

0.330 

0.364 

Midsection 

0.435 

256 

111.4 

0.80 

0.048 

0.026 

0.330 

0.404 

Outer 

0.230 

170 

39.1 

0.65 

0.026 

0.008 

0.120 

0.154 

Total 

- 

511 

238.4 

0.80 

0.090 

0.052 

0.780 

0.922 
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Fig. 1. - Electron bombardment ion thrustor. 
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Figure 2. - Electrode configuration showing the four groups of charge- 
exchange ions. 


Zero volt 
equipo- 
Screen tentiaU 



Neutraliza-I 
Primary ion tion plane-J 
beam boundary 


Screen 


(a) Typical trajectories. 



Aerator 

kV .^\ 


Neutraliza- 
tion plane-^ 
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(c) Energy distribution. 
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Figure 3. - Charge-exchange ion formation for group (2) ions with a 
central hole. Primary ion current of 1. 034 mA. Propellant utiliza- 
tion, 90 percent 



Figure 4. - Number of group (2) ions against dispersion angle. 










Figure 7. - Calculated values of neutral density for various dis- 
tances downstream of neutralization plane. 



Figure 8. - Calculated values of group (4) charge-exchange ions 
passing through the surface of a 15-centimeter diameter cyl- 
inder, the assumed primary ion beam boundary. 
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(b) Dispersion angle, deg 

Figure 9. - Two surfaces considered as intercepting thrus- Figure 10. - Charge-exchange ion flux 

tor fluxes. arriving at the two surfaces shown in 

figure 9. 



